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Electronic structure and magnetic properties for iron deficient TlFe2-ccSe2 compounds are studied by first- 
principles calculations. We find that for the case of x = 0.5 with a Fe-vacancy ordered orthorhombic super- 
structure, the ground state exhibits a stripe-like antiferromagnetic ordering and opens a sizable band gap if the 
short-ranged Coulomb interaction of Fe-3d electrons is moderately strong, manifesting a possible Mott insu- 
lating state. While increasing Fe-vacancies from the x = side, where the band structure is similar to that 
of a heavily electron-doped FeSe system, the Mott localization can be driven by kinetic energy reduction as 
evidenced by the band narrowing effect. Implications of this scenario in the recent experiments on TlFe2-ccSe2 
are discussed. 

PACS numbers: 71.10.Hf, 71.27.+a, 71.55.-i, 75.20.Hr 



The discovery of superconductivity fSC) with critical tem- 
peratures up to 56 K in iron pnictides luJ-la] has triggered re- 
newed interest in searching new route to high-temperature SC. 
Considerable concerns have been focused on the nature of the 
parent compounds which show variable bad metal behavior]^ 
and a universal strip-like antiferromagnetic (SDW) orderiS]. 
The magnetic ordering was proposed to be the consequence of 
the low energy states ( or itinerant electrons ) within the nearly 
nested Fermi surfaces, and hence the Fermi surface nesting is 
responsible to SC when the SDW order is suppressed 19-11]. 
An alternative possibility is that the magnetic structure is due 
to the strong correlation among the Fe-3d electrons, so that 
the states far away from the Fermi surfaces should be taken 
into account as well. To this end, Fermi surface nesting is not 
a necessary ingredient and the J1-J2 Heisenberg model based 
on the local moment picture could be an appropriate starting 

point. ilMilii] 

Most recently, Fang et al. (vk reported that the Fe-deficient 
compounds (Tl,K)Fe2-a:Se2 exhibit SC with Tc up to ~ 31 K 
forx = 0.12 ~ 0.3. Special interest in this class of materials 
is that the SC emerges in proximity to an insulating phase. It 
is yet unknown whether this insulating phase is a Mott insu- 
lator driven by Fe-3d electron correlations. Theoretically, this 
possibility is mysterious, as previous first-principle calcula- 
tions on both TlFesSes fls*] and KFesSes fl?, ^^with the 
ThCr2Si2 structure (122-type) suggest that the parent com- 
pounds of the ternary iron chalcogenides should be metallic 
with either checkerboard antiferromagnetic (AFM) (for Tl- 
122) [18] or SDW (for K- 122) [20] order, much like the elec- 
tron overdoped 11 -type iron selenides. Experimentally, it has 
been reported recently that the alkali intercalated compounds 
Ko.8Fe2Se2L21] and Cso.8Fe2Se1.96 [22] ( both iso- structural 
to BaFe2 AS2), while superconducting under 30 K and 27 K re- 
spectively, exhibit the metallic behavior in their normal states. 

Two closely related questions thus arise: What is the ground 
state of the "parent" compound of these superconducting 
ternary iron chalcogenides of 122-type structure, and how can 
a Mott-insulating phase develop and then diminish with in- 



creasing Fe-content or electron doping giving way to SC? 

In this paper, we suggest partial answers to the questions 
by first-principles study on TlFe2-a:Se2. We start from and 
pay special attention to the case with x = 0.5, which is 
stoichiometrically equivalent to Tl2Fe3Se4 but with the 122- 
type structure. Hence, Fe^+ is the nominal valence as in 
other iron pnictides/chalcogenides. Our calculation shows 
that the ground state of TlFei.5Se2 with Fe-vacancy ordered 
orthorhombic superstructure is an SDW phase which can open 
a sizable gap if a moderately strong electron correlation is im- 
posed. The quantitative change of the bandwidth of Fe-3d 
electrons further supports the Mott localization driven by ki- 
netic energy reduction due to Fe-vacancies. This scenario is 
quite similar to the Mott localization proposed for iron oxy- 
chalcogenides LaoQoFeoOfSe.S^o liill . where the reduction of 
kinetic energy or the band narrowing is due to the expanded 
interatomic Fe-Fe distance. We will discuss the implications 
of this scenario in the recent experiments on TlFe2-xSe2. 

The electronic structure calculations were performed with 
the Vienna Ab-initio Simulation Package (VASP)!!!!!!. All 
structures were optimized so that the forces on individual 
atoms were smaller than 0.02 eV/A and the pressure conver- 
gence criterion is chosen to be 0.5 kbar. For the optimization 
and ground state calculations, a 8 x 4 x 4 Monkhorst-Pack 
k-grid[26] was employed, while 16 x 8 x 8 Monkhorst-Pack 
k-grid was used for the density of states (DOS) calculations. 
The PBE flavor of general gradient approximation (GGA) to 
the exchange-correlation functional [27] was applied through- 
out the calculations. 

We first discuss the possible crystal structure and spin con- 
figurations of TlFei.5Se2. The Mossbauer experiment sug- 
gests a body-centered orthorhombic (BCO) crystal structure 
[28], where the Fe sheets form superstructure due to Fe- 
deficiency (FIG.[T]). There are four possible stacking config- 
urations for the two neighbouring Fe-layers[32], as shown in 
FIG.[T1 We consider the AFM and SDW configurations within 
each layers, while the interlayer coupling along c-axis can be 
either ferromagnetic (FM) or AFM. The relative energies for 
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FIG. 1: Superstructure of TlFei.5Se2. A, B, and C correspond to 
three different sites directly above the Fe- vacancy, the 3 -coordinated 
Fe site, and the 2-coordinated site, respectively. Another stacking 
pattern A' is that the upper-layer is rotated 90° around c-axis. 



TABLE II: Geometry and magnetic properties of TlFei.5Se2. Only 
the ground state (SDW^) configurations are listed, zse is the inter- 
nal coordinate for Se and mpe is the local magnetic moment on Fe 
atoms. The numbers outside and inside the brackets for mpe are for 
3-coordinated and 2-coordinated Fe atoms, respectively. 







AA 


AB 


AC 


SDW" 


a 


5.4824 


5.5407 


5.5483 




b 


11.0519 


11.0497 


11.0270 




c 


13.8162 


13.7892 


13.8176 




ZSe 


0.351 


0.351(0.350) 


0.350 




mFe(Ms) 


2.60(2.73) 


2.61(2.74) 


2.62(2.74 ) 



different stacking and ordering patterns are listed in Table Jl 
We find that the ground state ( for x = 0.5 ) should be of 
essentially two-dimensional SDW, as the energy differences 
among the different inter-layer magnetic orderings are very 
small (< 5meV/Fe), indicating the negligible inter-layer mag- 
netic coupling at this state. Hence, without losing generality, 
we shall focus on the AA- stacking pattern in the following 
discussions. 

We also notice that for the NM state, the full structural op- 
timization leads to the well-known c-collapse problem; while 
for the SDW state, the resulting structure (c and zse) is within 
the expectation. Thus, unless otherwise specified, the elec- 
tronic structure calculations in the following discussions were 
performed with the geometry relaxed in the SDW configura- 
tion. The calculated local magnetic moment of Fe atom is then 
mpe ~ 2.6^2.7 jj^B- 

TABLE I: Configuration energies of TlFei.5Se2 ( in unit of meV ). 
AA, AB, AC and AA' indicates the stacking geometry, while NM, 
AFM, SDW indicate the intralayer magnetic configurations. The su- 
perscripts a and / indicate interlay er AFM and FM orderings. 





AA 


AB 


AC 


AA' 


NM 


0.0 


9 


3 


12.7 


AFM-^ 


-200 


-111 


-172 




AFM" 


-173 


All 


-174 




SDW^ 


-268 


-268 


-270 




SDW" 


-273 


-268 


-269 





We now examine the DOS of the AA-stacking SDW^ plot- 
ted in Fig. |2(a)[ which exhibits a sharp dip around the 
Fermi energy. We find that the band gap Eg is vanish- 
ingly small for the fully optimized system under the GO A 
method. Furthermore, the AB- and AC-stacking SDW^ 
states, which are almost degenerate to the AA-stacking 
SDW^, are both metallic, with 0.216 states/(eV -Formula) 
and 0.656 states/(eV -Formula) respectively. Hence, the three 
stacking patterns may coexist at low temperatures. The SDW 
band gap, if exists, is too small to compare with the activation 
gap Ea ~ 57.7 meV observed in the x = 0.5 sample IitI] . It 
indicates that the observed activation gap in TlFe2-a:Se2 ( at 
least for X = 0.5 ) is not due to the SDW ordering itself. 
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FIG. 2: Total and projected density of states of TlFe2-ccSe2. |2(a)| 
The [/-dependent DOS at the SDW state for x = 0.5; [2(E)1 The 
x-dependent DOS at the NM state for [/ = with Fe-vacancy su- 
perstructures. In |2(a)| only a- spin DOS is shown since both spin are 
degenerate. All plots are renormalized to per TlFe2-a;Se2 formula. 
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Here, we suggest that the sizable activation gap can be 
attributed to the Mott localization driven by the moderately 
strong electron correlation. To seek for this possibility, we ex- 
tended our calculations by using the GGA+/7 method. The 
calculated DOS for = eV and U = 4.0 eV are compared 
in FIG. |2(a)( Within GGA+/7, an insulating gap develops im- 
mediately with increasing U, which turns out to be 40, 60, 80, 
140, and 230 meV for /7=1,2,3,4, and 5 eV, respectively. The 
projected DOS associated with Fe-3d electrons move to higher 
energies mainly distributed around 1 eV above the Fermi en- 
ergy ( for /7=4 eV ). 




vide an indication for the transition at certain < Xc < 0.5, 
i.e., when x > Xc, the Mott localization takes place. As 
the Fe-Fe distance increases less than 1% from x = to 
X = 0.5 ifivll . this amount of lattice expansion is not sufficient 
for the Mott localization in TlFe2-icSe2, as compared to iron 
oxychalcogenides La202Fe20(Se,S)2l23]. Here, we argue 
that the kinetic energy reduction caused by the Fe- vacancies 
should play an crucial role in driving the system to the insulat- 
ing phase. An intuitive estimate is that the coordinate number 
of Fe in the TlFei.5Se2 is reduced to 3 or 2 depending on the 
Fe site comparing to 4 in a perfect square lattice. Thus the to- 
tal kinetic energy is substantially reduced by the Fe- vacancies, 
enhancing the normalized electron correlation U/W, with W 
being the bandwidth proportional to the kinetic energy. 

Our argument can be checked by fitting band structures of 
X =0, 0.4 and 0.5 to the tight-binding models. We only need 
to consider the NM state ( and U = 0) and calculate the sum 
of the nearest neighbor hoppings ( absolute values ) among 
all five d-orbitals around a specific Fe atom. It serves as an 
approximate upper limit of the total kinetic energy Ek . Us- 
ing the fitted hopping parameters we obtain Ek =15.59 and 
9.98 eV for x =0 and 0.4, respectively. While for x =0.5, 
we obtain Ek =11.34 or 6.87 eV, corresponding to 3- or 2- 
coordinated Fe- sites, respectively. As the ratio of their num- 
bers is 2: 1, the average kinetic energy is 9.85 eV. For compari- 
son, the DOS at the NM state when [/ = for x = 0, 0.4, and 
0.5 are plotted in Fig. |2(b)[ The Fe-3d bandwidths are roughly 
4.8 eV, 3.8 eV, and 3.5 eV, respectively. Thus, we indeed find 
a substantial enhancement of the normalized electron correla- 
tion, manifesting the kinetic energy reduction caused by Fe- 
vacancies. 



FIG. 3: The Fermi surfaces ( for x = 0.5 and [/ = ) reconstructed 
using the MLWFs shown in the F-centered reciprocal BCO lattice. 



We also performed calculations for x = 0, 0.4, keeping the 
Fe-vacancy ordered tetragonal superstructure in the latter case 
lll7l|28|]. For x = the AFM state becomes more stable than 
the SDW state ( in agreement with Ref. llTsI] ), while the over- 
all band structure and Fermi surface ( FS ) of TlFe2As2 are 
similar to that of KFe2Se2 [20]. For illustration, the band 
structures for x = 0, 0.5 in the NM states are fitted to a tight- 
binding model Hamiltonian by using the maximally localized 
wannier function (MLWF)|29, 30] method. The results are 
plotted in Fig. (H where the color indicates the percentage 
of dzx{y) and dy2_^2 composition, from 0% (blue) to 100% 
(red). For x = 0, the Hamiltonian is particularly simple with 
5 -bands. The nearest neighbor hopping parameters and the 
on-site energies are listed in TAB.HVl The bands near Ef are 
dominated by dzx{y) and d^2_^2 orbitals except for one empty 
band around X. For x = 0.5, the Hamiltonian is too compli- 
cated due to the structural distortion, and all 5 d-orbitals are 
entangled considerably around Ep. Nevertheless, the FS of 
TlFei.5Se2 can be reconstructed as shown in Fig. O 

The apparent different band structures for x = and 0.5, 
together with the drastic change in their FS topologies, pro- 



TABLE III: Pressure dependency of crystal structure and band gap 
for X = 0.5 and U = 4.0 eV. 



P (GPa) 





2 


4 


6 


a (A) 


5.5507 


5.4517 


5.3716 


5.3079 


b(A) 


11.0342 


10.8842 


10.7283 


10.6177 


c(A) 


13.8027 


13.5166 


13.2874 


13.1102 


Eceii (eV) 


-187.01 


-186.65 


-185.91 


-158.00 


Eg (meV) 


140 


80 


50 






Finally, we remark that our calculations assume homoge- 
nous formation of the Fe-vacancy ordered superstructures. If 
this is the case in real materials, the activation gap observed 
in the X = 0.5 sample should be not due to the Anderson 
localization caused by vacancy disorder effect. As the Mott 
(de-)localization is sensitive to physical pressure which can 
lead to monotonous bandwidth expansion by lattice contrac- 
tion, we expect that the gap dependence on pressure can be 
used to test this scenario. The numerical results for the pres- 
sure dependent band gap Eg with fixed U = 4.0 eV are listed 
in TAB.HIIlf X = 0.5). From these results we expect that the 
activation gap, which is about 60 meV or well below 150 meV 
1 17] under ambient pressure, can be completely suppressed by 
applying pressure up to ~6 GPa and then the SC may emerge. 
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FIG. 4: Band structures of TlFe2-xSe2 fitted using MLWF. 



TABLE IV: Tight-binding Hamiltonian for TlFe2Se2. Only the near- 
est neighbor hoppings and on-site energies ( diagonal line in the 
brackets ) are shown here. All numbers are in eV. 





d^2 


dzx 


dzy 


d^2_y2 


dxy 




-0.019 
(4.765) 





0.156 


0.303 





dzx 





-0.028 
(5.084) 








0.286 


dzy 


0.156 





-0.314 
(5.084) 


-0.344 







-0.303 





-0.344 


-0.378 
(4.600) 


0.0 


dxy 





0.286 








-0.047 
(5.056) 



To summarize, we find that for TlFei.5Se2, the electronic 
band structure shows a sizable band gap if the short-ranged 
Coulomb interaction U beyond LDA/GGA is considered. 
From the experimentally observed activation gap, U should be 
at least 2eV. The corresponding ground state is then a stripe- 
like anti-ferromagnetic Mott insulator. The superstructure of 
various Fe- vacancy ordering patterns is important for the sta- 
bility of the Mott insulating phase and kinetic energy reduc- 
tion caused by the vacancies plays a crucial role for the Mott 
localization. The pressure dependence of the gap behavior 
could be used to test this possibility. 
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Note added: After completing this work, we became aware 
of two recent papers by X.W. Yan et al.[31], where indepen- 
dent first-principles studies on the cases of x = and x = 0.5 
were reported. 
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